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Abstract 

We perform in this paper a statefinder diagnostic to a dark energy model with two scalar fields, 
called "quintom" , where one of the scalar fields has a canonical kinetic energy term and the other 
has a negative one. Several kinds of potentials are discussed. Our results show that the statefinder 
diagnostic can differentiate quintom model with other dark energy models. 
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It is widely believed nowadays that the present universe is undergoing an accelerating 
expansion. The converging evidences come from the analysis of data from supernova 
CMB S, S fl J and WMAP Q. In order to explain the cosmic accelerating 

expansion, a modified theory of gravity or the existence of dark energy is required. Perhaps 
the simplest candidate of dark energy is the cosmological constant with equation of state w = 
p/p = —1. However, there exist two problems with it, i.e., why the cosmological constant is 
so tiny compared to the theoretical expectations and why the remnant cosmological constant 
is becoming visible precisely at the present time in the cosmic history and why it does not 
exactly vanish. The inspiration coming from inflation has suggested that dark energy models 
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. The quintessence scalar field has a positive kinetic 
term in its Lagrangian, which violates the strong energy condition but not the dominant 
energy condition, and the evolution of its equation of state parameter w is in the range of 
— 1 < w < 1. The phantom scalar field on the other hand possesses a negative kinetic term in 
its Lagrangian, which leads to some strange properties, such as the violation of the dominant 
energy condition and the occurrence of Big rip [66J, |67|, |§a |6Ja, |70j, LLU IZ-J l22l IZfa |75J] ( Let 
us note here that the possible avoidance of the Big Rip by various kinds of physical effects 

has been discussed by many authors 0, 00 0, Q)- 

The evolution of 

the equation of state parameter w for phantom is in the rang e of w < —1. Other models 



including the Chaplygin gas 



models 



00 



brane world models 



H 0, 3, 



holographic 



9l| et al are also proposed to account for the present accelerating cosmic expansion. 



The recent analysis of various cosmological data seems to indicate that it is mildly favored 
that the evolution of the equation of state parameter of the dark energy changes from 



w > — 1 to «; < - lat small redshift 921. 



state parai 
0,HQ. 



Although this result may be considered 



as tentative only since the standard A CDM model remains inside ~ 2a error bars for 
all data, it is worthwhile to explore theoretical models which can bring w crossing- 1. In 
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this regard, most models mentioned above can not do the job. "Quintom" 
which assumes that the dark energy is composed of quintessence and phantom, on the other 
hand, can implement w crossing -1 and in some cases fits the cosmological data better than 



models with w > —1. In addition "quintom" model can reconcile the coincidence problem 
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and predicts some interesting features in the evolution and fate of our universe 
Considering a flat homogeneous and isotropic FRW universe filled with matter, radiation 
and "quintom" dark energy which is composed of a phantom scalar field and a normal 
scalar field <p, the dynamic equations for the universe, phantom and normal scalar fields can 
be described by 

^matO ^radO 
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(1) 



+ 3-0-^(0,^=0, (2) 

(X 

<p + 3% + V (<f>, tp), (3) 

where V(0, (p) is the interacting potential between phantom and quintessence, denotes the 
present time, a dot means a derivative with respect to time t and a is the scale factor. The 
energy density and pressure for the quintom can be expressed as 

p q m = -0 2 + 2 + V(0,y?) , (4) 

Pqui = ~0 2 + V 2 - V{<f>, V) • (5) 

Thus we can obtain the equation of state for the quintom 

Pqui -0 2 + 2 - V{(f), (f) , , 

W au i = — — = = • 6 

Pqm -0 2 + 2 + n<^) 

This expression implies w qui > — 1 when \<p\ > |0| and w qui < — 1 when \<p\ < |0|. 

Since there are more and more models proposed to explain the cosmic acceleration, it 
is very desirable to find a way to discriminate betw een t he various contenders in a model 

nni 

independent manner. In this regard, Sahni et al. [99, 100] recently proposed a cosmo logical 
diagnostic pair {r, s} called statefinder, which are defined as 

r = ^IP' S= 3( g -l/2) ' [ V 

to differentiate between different forms of dark energy. Here q is the deceleration parameter. 
Apparently statefinder parameters only depend on a, and is thus a geometrical diagnostic. It 
is easy to see that statefinder differentiates the expansion dynamics with higher derivatives 



of scale factor a and is a natural next step beyond H and q. Since different cosmo logical 
models involving dark energy exhibit qualitatively different evolution trajectories in the s — r 
plane, this statefmder diagnostic can differentiate various kinds of dark energy models. For 
LCDM (or ACDM) cosmological model, which consists of a mixture of vacuum energy and 
cold dark mass, the statefmder parameters correspond to a fixed point {r = 1, s = 0}. By 
far some models, including the cosmological constant, quintessence, phantom, the Chaplygin 
gas ; br aneworld mo dels , holog r aphi c models, interacting and coupling dark energy models 
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1061 ] . have been successfully differentiated. For example, 
although the quintessence model with inverse power law potential, the phantom model 
with power law potential and the Chaplygin gas models all tend to approach the LCDM 
fixed point, for quintessence and phantom models the trajectories lie in the regions s > 0, 
r < 1 while for Chaplygin gas models the trajectories lie in the regions s < 0, r > 1, 
and on the other hand the quintessence tracker models and the Chaplygin gas models have 
typical trajectories similar to arcs of a parabola (u pward and downward ) respectively while 
for phantom model the trajectories are different 0, Eml, Ed3, 
quintessence models the trajectories of r(s) form swirl before reaching the attractor [105 1. 

In this paper we apply the statefmder diagnostic to the quintom dark energy model. To 
begin with, let us use another form of statefmder parameters which can be written as 

9{p + p)p {p + p)P fa . 

2 P P P P 

Here p is the total energy density and p is the total pressure in the universe. Since the total 
energy, quintom energy and radiation are conserved respectively, we have p = —3H(p + p), 
Pqui = -3H(1 + w qui )p qui and p rad = -AHp rad . Thus we can obtain 

3 

r = 1 - -[wgui/H - 3w qui (l + w qui )]VL qui + 2Q rad , (9) 
g _ -3[w qui /H - 3w qui (l + w qui )}Q qui + 4Q rad 

and 

q = ~(1 + 3w quj Vt qui + Vt rad ) . (11) 
Here Vt qui = p qui /p and Vt rad = p rad /p. 
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In the following we will discuss the statefinder for the quintom model with several kinds 
of potentials. Firstly we assume that there is no direct coupling between the phantom scalar 
field and the normal scalar field with such a potential V(<f>, <p) = V^oe~ a ^ + V^eT^ ', where 
a and (3 are constants. In this case the universe has the phantom dominated late time big 
rip attractor j^]. In Fig. ^ we show the time evolution of statefinder pair {r, s} in the 
time interval — G [0.5,4] where to is the present time. The model parameters are chosen as 
V^o = 0.3p and V^ = 0.6p , where p is the present energy density of our universe. We 
see that in the past and future the r — s is almost linear, which means that the deceleration 
parameter changes from one constant to another nearly with the increasing of time, and the 



parameters will pass the fixed point of LCDM in the future. 



different from other dark energy models discussed in Refs. 
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FIG. 1: The r — s diagram of no direct coupling exponential potential V(4>, ip) = V^e - "^ + V ip oe~ /3lfi . 
Curves r(s) evolve in the time interval j- G [0.5,4] where to is the present time. In the left figure 
P = 3 and a = 1, 1.5, 2(solid line, dashed line and dot-dashed line respectively). In the right figure 
a = 1 and (3 = 2, 2.5, 3(solid line, dashed line and dot-dashed line respectively). Dots locate the 
current values of the statefinder parameters. 



Next, two cases of potential V(4>, ip) = V^e - ^ + V^e'^ + V e- K ^ + ^ and V(4>, ip) = 
V^oe -Q ^ 2 + V^oe - ^ 2 will be studied, as before where a, (3 and k are constants. The quintom 
model with both above potentials have the late time attractors but t he fo rmer potential 
leads to a big rip attractor and the latter to a de Sitter attractor 

• In Fig. H the 

curves of r(s) are given. Apparently the left figure is very similar to Fig. ^ This shows that 
for uncoupling and coupling exponential potentials the evolutions of our universe are very 
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similar in the time interval we consider here. The right figure is differe nt fr om Fig. ^ but 
has a common characteristic with the phantom with power la w po tential jl03^ , quintessence 
with inverse power law potential and Chaplygin gas model 10l| that it reaches the point 
of LCDM with the increasing of time. The reason is that they all lead to the same fate of 
the universe-de Sitter expansion, but the trajectories to LCDM are different, therefore they 
can be differentiated. 
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FIG. 2: In the left figure curves r(s) evolve in the time interval j- G [0.5,4] where to is the present 
time for the potential V(<j), ip) = V^e'^ + V {fi0 e~ (3ip + V e- K( - 4>+ip \ The model parameter are chosen 
as V^o = 0.3po) VipO = 0.6po> Vq = 0.3/9o, a = 1, [i = 1 and k = 1, 2(solid line and dot-dashed 
line respectively). In the right figure the potential is V(<p, ip) = V^eT^ + V^e - ^ 2 . The model 
parameter are chosen as V^o = 0.3po, V^q = 0-Qpo, a = 1 and = 1, 3(solid line and dot-dashed 
line respectively). Dots locate the current values of the statefinder parameters. 



Finally, we turn to the case of linear coupling potential V((f), if) = + ip) + A0y9, where 
k an d A a re two constants. The scalar field with a linear potential was firstly studied in 
Ref. 108^ and it has be e n ar gued that such a potential is favored by anthropic principle 



considerations 



109 



110 



1 1211 . In addition if 



and can solve the coincidence problem 
the universe is d omin ated by quintessence (phantom) with this potential it ends with a big 
crunch(big rip) 113j . The time evolutions of statefinder pair {r, s} in the time interval 
Y G [0.5,4] are shown in Fig. El We see that in the case of negative coupling the diagram 
is very similar to Fig. ^ and the left figure in Fig. which shows that in these cases the 
evolutions of our universe are similar in the time interval we consider here, as shown in 
Fig. |3J But for A > the figure of {r, s} is very different. In order to well explain this we 
give the time evolutions of {r, q} in Fig. From Fig. 0] and Fig. El we find if the coupling 
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constant A is positive the accelerating expansion of our universe is temporal and the universe 
will decelerate again. 




s s 

FIG. 3: The r — s diagrams of linear potential case. Curves r(s) evolve in the time interval 
€ [0.5,4]. In the left figure k = 0.8 and A = —0.2, — 0.8(solid line and dashed line respectively). 
In the right figure k = 0.8 and A = 0.2, 0.8(solid line and dashed line respectively). Dots locate 
the current values of the statefinder parameters. 
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FIG. 4: The evolution of scale factor a with t/to. 1 an d 2 represent potentials V(4>, if) = V^oe + 
V^e-^ and V(<f>,<p) = V^e"^ + V v0 e~^ + Vbe~ K ^+^ respectively. 3 and dot dashed line denote 
the linear potential with negative and positive coupling respectively. 

In summary, we study the statefinder parameters of the quintom dark energy model. 
Several kinds of potentials are discussed. It is found that the statefinder diagnostic can 
differentiate the quintom model with other dark energy models, but seems not very helpful 
to differentiate quintom dark energy models with some different kinds of potentials which 
lead to a similar evolution of our unverse in the time interval we consider. 
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FIG. 5: The r — q diagram of linear potential case. Curves r(q) evolve in the time interval 
€ [0.2,4], k = 0.8 and A = 0.2, 0.8(solid line and dashed line respectively). Dots locate the 
current values of the statefinder parameters. 
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